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Abstract 

In this paper we show an 0{in?^^ log^ n) time algorithm for finding a maximum flow 
in a planar graph with multiple sources and multiple sinks. This is the fastest algorithm 
whose running time depends only on the number of vertices in the graph. For general (non- 
planar) graphs the multiple-source multiple-sink version of the maximum flow problem is as 
difHcult as the standard single-source single-sink version. However, the standard reduction 
does not preserve the planarity of the graph, and it is not known how to generalize existing 
maximum flow algorithms for planar graphs to the multiple-source multiple-sink maximum 
flow problem. 

1 Introduction 

In the standard maximum flow problem, we wish to send the maximum possible amount of flow 
from a specific vertex designated as a source to another vertex designated as a sink^ this is the 
single-source single-sink maximum flow problem. In a more general formulation of the problem 
we have many sources and many sinks, this is the multiple- source multiple-sink maximum flow 
problem. Ford and Fulkerson suggested a simple reduction from the multiple-source multiple- 
sink maximum flow problem to the standard single-source single-sink problem. We add to the 
graph a new vertex s' (called super-source) and connect it to all the sources in the graph with 
arcs of infinite capacity. Similarly we add a new vertex t' {super-sink) and connect all the 
sinks in the graph to it with arcs of infinite capacity. Then we solve a single-source single-sink 
maximum flow problem with s' as the source and t' as the sink. 

In this work, we are interested in flow networks whose underlying graph is planar. The 
fastest known algorithm for the single-source single-sink maximum flow in planar graphs is the 
algorithm of Borradaile and Klein [T], which runs in O(nlogn) time, where n is the number 
of vertices in the graph. This algorithm was simplified by Schmidt et al. [19J and by Erickson 
See [1] for a survey on previous results for the problem. However, it is not known how 
to generalize existing planar maximum flow algorithms to the multiple-source multiple-sink 
problem. The Ford and Fulkerson standard reduction does not preserve the planarity of the 
graph. Consider a planar graph in which four sources form a clique in the underlying undirected 
graph. When we add the super-sink s' to the graph we get a complete graph of 5 vertices in the 
underlying undirected graph, which is not planar. This standard reduction, using an algorithm 
for maximum flow in general graphs such as the algorithm of Goldberg and Tarjan [6], gives us 
an 0(n^ log n) time algorithm for our problem. 

In addition to uses in standard applications of planar graph, such as transportation net- 
works, the multiple-source multiple-sink maximum flow problem in planar graphs has many 
applications of energy functions minimization in computer vision. The maximum flow problem 
gives us a solution to the related minimum s — t cut problem [5j. Greig et al. ^ were the first 
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to use minimum s — t cut algorithms for such a minimization problem of an energy function, 
for binary image restoration. Since then, many other similar uses have emerged in computer 
vision, such as image segmentation, multi-camera stereo vision, object recognition, and others, 
see [3j for a survey. From a graph point of view, the graphs for these applications are all similar. 
For every pixel in the image, the graph contains a vertex, such that vertices of two neighboring 
pixels are connected by an edge. In addition, we add a source s and a sink t and connect them 
to the vertices of the pixels. The pixels define a planar graph, in fact a grid, but the source and 
the sink destroy the planarity of the graph. However, by breaking s into multiple sources and 
t into multiple sinks, we can reduce the graph to a multiple-source multiple-sink planar graph. 

Miller and Naor [IT] studied the problem of maximum flow with multiple sources and multi- 
ple sinks in planar graphs. When the demand (the required difference between the incoming and 
outgoing amounts of flow) at each vertex is known, they showed how to reduce this problem to 
a circulation problem, which can be solved in planar graphs in 0(nlog^n/loglogn) time [18] . 
Miller and Naor also showed an 0(n log^'^^ n) time algorithm for the case where all the sinks and 
the sources are on the boundary of a single face, and generalized it to an 0(A:^n^/^ log^ n) time 
algorithm for the case where the sources and the sinks reside on the boundaries of k different 
faces. The time bound of the first algorithm can be improved to O(nlogn) using the linear time 
shortest path algorithm of Henzinger et al. [10], and the time bound of the second algorithm 
can be improved to 0(A;^nlog^n) using the 0(n log n) time single-source single-sink maximum 
flow algorithm of Borradaile and Klein ^ . An additional approach that Miller and Naor show 
iterate over all the sources and all the sinks, and solve the maximum flow problem from every 
source to every sink, each time taking the residual network of the previous flow as an input. 
This gives an 0(|5'||T|nlogn) time bound, where 151 is the number of sources and \T\ is the 
number of sinks, using the algorithm of Borradaile and Klein [T|. 

Recently, Borradaile and Wulff-Nilsen [2j and Klein and Mozes [13] showed two algorithms for 
maximum flow in planar graphs with multiple sources and a single sink (or equivalently, a single 
source and multiple sinks) in 0(n^/^ logn) time. This leads to an ©(minllS"!, |T|}n'^/^ logn) time 
algorithm for our problem. 

In this paper, we show an 0{n^^'^ log^ n) time algorithm for the multiple-source multiple-sink 
maximum flow problem in planar graphs. This is the fastest algorithm for the problem whose 
running time depends only on the value of n, and not on the set of sources, the set of sinks, or 
the capacities of the arcs. 

Both algorithms for the multiple-source single-sink version of the problem [21 [T3] use pre- 
flows to find the required flow, we use the similar concept of pseudoflows for our algorithm. 
Pseudoflows were originally used for the minimum cost flow problem [7J, Hochbaum [llj was 
the first to use pseudoflows for the maximum flow problem. As the two algorithms of [3 I14j . 
we also decompose the problem into smaller subproblems using cycle separators [16]. Cycle 
separators were first used for maximum flow in planar graphs by Johnson and Venkatesan [TS] . 
We use an approach similar to the one of Miller and Naor [Ij] in our algorithm, initially our 
flow function / is all zero, then we repeatedly solve some maximum flow problem in the residual 
network with respect to / and add to / the resulting flow. 

2 Preliminaries 

We consider a simple directed planar graph G = {V,E), where V is the set of vertices and E 
is the set of arcs. We denote the number of vertices by n, since the graph is planar we have 
\E\ = 0{n). For v G V, in{v) = {{u,v) \ {u,v) G E} is the set of incoming arcs. We assume 
that G is given with a fixed planar embedding, in other words it is a plane graph. Combinatorial 
representation of such an embedding can be found in 0(n) time [12j. We assume that the graph 
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is connected, as we can process every connected component separately. 

A flow network consists of a graph G, a set S of vertices designated as sources, a set T of 
vertices designated as sinks, and a capacity function that assigns to every arc e € E a finite 
capacity c(e) > 0. A function f : E ^ Ris a flow function if and only if it satisfies the following 
three constraints: 

ye £ E (capacity constraint), (1) 

\/{u,v) G E (antisymmetry constraint), (2) 

\/v G V \ S U T (flow conservation constraint). (3) 

For the antisymmetry constraint we assume that for every arc {u,v) in E, the arc {v,u) is also 
in E, if this is not the case to begin with, then we add the arc with capacity 0. In the planar 
embedding of G we identify the arcs {u,v) and {v,u), that is, we embed them as a single edge. 
In other words, an edge is a pair of antiparallel arcs. 

The value of a flow / is 'Yle&m{t)t<^T fi^)-: ^he amount of flow that enters the sinks. A 
maximum flow is a flow of maximum value. 

Goldberg and Trajan [7J define pseudoflow by removing the flow conservation constraint, in 
other words, a pseudoflow is a function f : E ^ R that satisfies the capacity constraint and 
the antisymmetry constraint. The flow excess of a vertex v is excess{v) = Yleein{v)) /(^)- 
can view any pseudoflow / as a flow function, if we choose the right sets of sources and sinks. 
Indeed, if we replace S with S' = {v \ excess{v) > 0} and T with T' = {v \ excess (v) < 0} then 
/ is a flow from S' to T' with respect to the graph G and the capacity function c. 

The residual capacity of an arc e with respect to a pseudoflow / is rj(e) = c(e) — /(e). An 
arc e is residual if rj(e) > 0. A residual path is a path of residual arcs. A flow is maximum 
if and only if it has no residual path from 5 to T [5]. The residual network with respect to 
a pseudoflow / is the flow network G/ on the graph G with sources S, sinks T, and capacity 
function rj. 

2.1 Cycle separators 

Assume we assign weights that sum to 1 to the vertices of G. Let C be a cycle in the plane 
(not a cycle of G) that meets G only in its vertices, and let P be the set of vertices in which G 
and G meet. The cycle C separates G into two edge-disjoint pieces - the subgraph of G inside 
C including the vertices of P, and the subgraph of G outside G including the vertices of P. If 
the total weight of vertices strictly inside C, and the total weight of vertices strictly outside C 
(we do not sum the vertices of P in the total) are both at most 2/3, then the set of vertices P 
is a cycle separator. 

Miller [16| showed that for any triangulated planar graph there is a cycle separator of 0{^/n) 
vertice^. Furthermore, such a separator can be found in 0(n) time. We can assume without 
loss of generality that G is triangulated, since otherwise we can triangulate it in linear time 
with arcs of capacity 0. 

For convenience, we use the term separator to describe a set of vertices that separates 
between two pieces even if the pieces do not satisfy the weight constraint. 

^Miller [IS] considers graphs that are not 2-connected as a special case which might have a single vertex 
separator instead of a cycle separator, for our purpose it does not matter. 
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3 Flow summation 

A key technique in our algorithm is flow summation, defined as fohows. Let / be a pseudoflow 
with respect to some capacity function c, and let /' be a pseudoflow with respect to the residual 
capacity rf. We let (/ + f'){e) = /(e) + /'(e). The function / + /' is a pseudoflow satisfying 
the capacity constraint with respect to c. Furthermore, if / and /' are both flow functions 
(satisfying flow conservation) from the same set of sources to the same set of sinks, then / + /' 
is also a flow function from this set of sources to the set of sinks. 

We use two algorithms of Miller and Naor [T7] to obtain a maximum flow through flow 
summation. In the first algorithm, we apply a single-source single-sink maximum flow algorithm 
for each pair of source and sink: 

Algorithm 1 

1. Set /(e) = for every arc e. 

2. For every s £ S do: 

2.1. For every t £ T do: 

2.1.1. Compute a maximum flow /' from s to t in the residual network Gf. 

2.1.2. Let /^/ + /'. 

3. Return /. 



Algorithm [T] is described in Section 7 of [T7] and also in Theorem 4.1 of |2j. The correctness 
of this algorithm implies the following lemma: 

Lemma 3.1. For a subset of the sources S' C S, let f be a maximum flow from S' to T, and 
let f he a maximum flow from S\S' to T in Gf. Then, f + f is a maximum flow from S to 
T. Symmetrically, for a subset of the sinks T' QT with a maximum flow f from S to T' and a 
maximum flow f in Gf from S to T \T' , the sum f + f is a maximum flow from S to T . 

The second algorithm of Miller and Naor |17j that we use is a recursive algorithm. For 
X C y we denote X n 5 by Xg, and similarly Xt is X r\T. The algorithm divides the vertices 
of the graph that are sources or sinks into two sets, L and R. After finding a maximum flow 
/ from Ls to i?t, the algorithm uses the cut with respect to / to decompose the problem into 
smaller problems. The cut with respect to / contains an arc {u, v) if and only if there is a 
residual path with respect to / form a vertex of Lg to u, but there is no residual path from Lg 
to V. After solving the maximum flow problem in the residual network inside each connected 
component, and adding the resulting flows to /, the algorithm finishes by computing maximum 
flow from Rg to Lt in the residual network. 

Algorithm 2 

1. Partition SUT into two disjoint sets L and R. 

2. Compute a maximum flow / from Lg to Rt, and let C be the cut with respect to /. 

3. Remove the edges of C from the graph. 

4. For each connected component K do: 

4.1. Find a maximum flow /' from Kg to Kt in the residual network Gf restricted to K. 

4.2. Let / ^ / + /'. 

5. Restore the edges of C to the graph 

6. Find a maximum flow /' from Rg to Lt in the residual network Gf. 

7. Let / ^ / + /'. 

8. Return /. 
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(a) Q' is a path carrying flow (b) Q' is a cycle of flow, 

from s G 5 to t e T. 

Figure 1: Proof of Lemma |3.2[ The path Q (dashed) is a residual path from a £ A to b £ B, 
which was created after adding Q' (solid). In both cases, when S A or T <^ B, we get that 
there was a residual path from A to B also before adding Q'. 



The correctness of Algorithm [2] is proven in Lemma 6.1 of [17j . It is described there for a 
specific choice of L and R, but as Miller and Naor note in Section 7 of [T7], the proof does not 
depend on this choice of the two sets. Intuitively, we can compute the maximum flow in each 
connected component separately because even if we had not removed the edges of C from G, 
then a residual path from a vertex in Kg to a vertex in Kt could not contain an arc of C, since 
these arcs are not residual. 

The following lemma allows us to use flow summation while keeping some invariant about 
non-existence of specific residual paths. 

Lemma 3.2. Let f be a pseudoflow such that there is no residual path with respect to f from 
any vertex in a set A <^V to any vertex in a set B CV , let f be a flow in Gf from S to T. If 
S ^ A or T (1 B, then there is no residual path from A to B with respect to f + f ■ 

Proof. We can view the flow /' as the sum of paths that carry flow from S to T and of cycles 
that carry flow [5]. Thus, we can view the summation of / and /' as adding the paths and 
cycles of flow that compose /' to / one by one in some arbitrary order. 

Assume for contradiction that the claim is not true. Before adding /' to / there was no 
residual path from A to B, then at some point a residual path Q was created from a vertex of 
A to a vertex of B. We choose Q to be the first such residual path created, as we added one by 
one the paths and cycles which compose /', and let Q' be the path or cycle of /' whose addition 
made Q residual (if the addition of Q' created more than one residual path from A to B, we 
choose one of them arbitrarily to be Q). 



Consider the case where Q' is a path (see Fig. 1(a)). The paths Q and Q' must share a 
common vertex v. If S" C j4 then we choose v such that v is the last vertex of Q' in Q, and let 
Q" be the path that begins with the prefix of Q' before v and ends with the suffix of Q that 
starts at V. If T C i? then we choose v such that v is the first vertex of Q' in Q, and let Q" be 
the path that begins with the prefix of Q before v and ends with the suffix of Q' that starts at 
V. In both cases, the path Q" was a residual path from A to B before we added Q' to the flow, 
contradicting the choice of Q. 



In the case where Q' is a cycle (see Fig. 1(b) ), let u he the first vertex of Q that is also a 



vertex of Q' and let v be the last vertex of Q common with Q' . Since Q' is a cycle it must 
contain a path from u to v. Let Q" he the path that we obtain by replacing the subpath of Q 
from n to f with the subpath of Q' . Again, the path Q" was a residual path from A to B before 
we added Q' to the flow, contradicting the choice of Q. □ 
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4 The algorithm 



The general framework of our algorithm is Algorithm [21 which we described in Section [3l This 
algorithm is recursive, for the base of the recursion we can use the cases where IS*! = 0(1), 
|T| = 0(1) or I^IITI = 0{y/n), in all of these cases we can find a maximum flow from S" to T in 
0(?7-^/^ log n) time, as described in the introduction. 

We use a procedure that computes a maximum flow between sources that are in one piece 
of the graph and sinks that are in another piece, with respect to a separator of size r. In the 
next section, we give an algorithm for this task which runs in Oir?!"^ logn + rn) time. 

We assign to each source or sink in G a weight of (l^l + |T|)~^, the other vertices have weight 
0, and find a separator P that separates the graph into two pieces - X and y, such that the 
total number of sources and sinks in each of them at most (2/3) (| 5*1 + |r|) + 0(^\/n) (the term 
0(-y/rz) comes from the fact that the vertices of P, that are common to X and y, can also be 
sources or sinks). We split the sources and sinks of G into two disjoint sets L and R (Step 1 of 
Algorithm [2]) as follows - L = Xg U \ Pg, i? = U Yj \ Pj . In other words, we define L and 
R according to X and y, and for the vertices of P, which are common to X and y, we put the 
sinks in L and the sources in R. 

Next (Step 2), we find a maximum flow from Lg to Rt- With the algorithm of Section [5l 
this takes 0(n^/^ logn) time. After we found the maximum flow /, it is easy to find its cut C 
in linear time. We remove the edges of C from the graph (Step 3). 

Now (Step 4), we find a maximum fiow inside each connected component. Let i^T be a 
component of G \ C with uk vertices. The component K is of one of two types, either K D Rt 
is empty or n is empty. 

For the first type we use Lemma [3. II to find a maximum fiow from Kg to Kt in three stages. 
First we find a maximum fiow from KCiRg (if it is not empty) to Kt = KOLf. Now it remains to 
find a maximum flow from KCiLg to Kt in the residual network. We split Kt into two parts, we 
find a maximum flow from KCiLg to KtCiPt in the residual network, and from KCiLg to Kt\Pt 
in the residual network (each time we compute a fiow, we add it to /). We use the algorithm of 
Section [5]to perform the first two maximum flow computations in 0{n^ logriK + |Pn K\nK) 
time (note that PDK remains a cycle separator between the sources and the sinks in these two 
cases), and we use a recursive application of our main algorithm for the third computation. For 
the second type of components we use Lemma l3.1l is a similar way, we compute maximum flow 
from Ks = KORg to KtCiLt, then from KgOPs to KtORt in the residual network, and last from 
Ks\Ps to KtCiRt in the residual network. Again, the first two maximum fiow computations are 
applications of Section O and the third is a recursive call. The different connected components 
are pairwise vertex-disjoint, so the total time in Step 4 for all the applications of the algorithm 
of Section[5]is 0(n^/^ logn). The sources and sinks in each recursive call are either all contained 
in L or all contained in R, and none of them is in P, so there are at most (2/3)(|S'| + |r|) sources 
and sinks in each call. 

We restore the edges of C (Step 5), and find a maximum flow from Rg to Lt in the residual 
network (Step 6), again this takes 0{n^^'^logn) time using the algorithm of Section [5l This 
concludes our implementation of Algorithm [2j 

The correctness of our algorithm is derived from the correctness of Algorithm [2] [T7j . We call 
the algorithm recursively on each connected component of G \ C, the components are disjoint, 
and the number of sources and sinks in each component is at most (2/3)(|S'| + |r|), therefore 
the total time bound of our algorithm is 0(n^/^ log^ n). 

Theorem 4.1. There is an 0{'n?/'^ log^ n) time algorithm that solves the multiple- source multiple- 
sink maximum flow problem in a planar graph with n vertices. 
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5 Maximum flow from one side of a separator to another 



In this section we present a solution for a special case of the multiple-source multiple-sink 
maximum flow problem, where all the sources are in one piece of the graph, and all the sinks 
are in another piece. We use this as a procedure in our main algorithm. This procedure is 
similar to a procedure used by Borradaile and Wulff-Nilsen [2] , but here we can use a multiple- 
source single-sink maximum flow algorithm (t2j or [14j ) as a sub-procedure, to avoid (explicit) 
recursion. 

Let X be the piece of G that contains 5, and let Y be the piece of G that contains T. Let 
P = pi, . . . ,pr be the separator of r vertices that separates between X and Y. Since P is a 
cycle separator, for every 1 < i < r it is possible to add an edge between pi and pi^i without 
violating the planarity of the graph, we will do so below. 

For simplicity we assume that no vertex of P is a source or a sink. If this is not true then we 
can fix it easily - assume that s G P is a source, then we add a vertex s' inside a face adjacent 
to s on the side of X of the separator, connect s' to s with an edge of capacity Ai, where M 
is the sum of all capacities of all arcs in G, and make s' a source instead of s. The case where 
t G P is a sink is handled symmetrically. Note that this transformation keeps the number of 
vertices and the number of edges in the graph 0{n). 

We begin by finding a maximum flow fx from S to all vertices of P, inside the piece X. 
We add to X a super-sink t in the place where Y is in the plane, and connect all the vertices 
of P to t with arcs of capacity A4. Since P is a cycle separator, adding t to X preserves the 
planarity of X. We find a multiple-source single-sink maximum flow from S to t, then fx is the 
restriction of this flow to X. Symmetrically, we find a maximum flow fy from all the vertices 
of P to T, inside the piece Y. 

Let / = fx + fy- Since X and Y are edge disjoint, / is a pseudoflow in G (that is, it satisfies 
the capacity constraint and the antisymmetry constraint). In the pseudofiow / the vertices that 
are not sources and not sinks that may have non-zero flow excess are only vertices of P. 

Next we show how to transform / from pseudoflow to a flow. We do so by balancing the 
flow excess of each vertex of P with non-zero excess, one vertex after the other, while keeping 
the following invariant: 

Invariant 5.1. Let P' <^ P be the vertices of P that we did not process yet. There is no residual 
path with respect to f from a vertex of S to a vertex of T, or from a vertex of S to a vertex of 
P' , or from a vertex of P' to a vertex ofT. 

Initially the invariant is true, there are no residual paths from 5* to P because fx is maxi- 
mum, and no residual paths from P to T because fy is maximum, also any path from S to T 
must contain a vertex of P, so there are no residual paths from S to T. When we are done, P' 
is empty, and all the vertices of P have flow excess, so / is a flow function. Since there is no 
residual path with respect to / from 5 to T, / is a maximum flow. 

We go over the vertices of P, from p\ to pr- Let pi be the current vertex that we process, 
P' = {Pj \ j ^ i} is the set of vertices that we did not process yet. If excess{pi) = we do 
nothing and proceed to the next vertex. 

The second case is when excess (pi) > 0. In this case, if i < r, we add arcs {pr,Pr-i), 
{pr-i,Pr-2), • • • ) iPi+2,Pi+i), each with capacity Ai. We flnd a maximum flow /j, whose value 
is bounded by excess{pi), from pi to pi^i in the residual network. We bound the value of the 
maximum flow by excess {pi) by adding a source s' inside the face that is adjacent both to pi 
and to Pi+i and connecting s' to pi with an arc of capacity excess{pi), then we actually flnd a 
maximum flow from s' to Pi+i- We remove the arcs we added from pr to Pi+i from the graph, 
and let / = / -|- /j. This takes 0(n) time using the algorithm of Hassin [9] for maximum flow 
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in a planar graph where the single source and the single sink are on the same face, with the 
shortest path algorithm of Henzinger et al. [10]. 

Lemma 5.2. After adding fi, the pseudoflow f satisfies Invariant \5.1l 

Proof. As in the proof of Lemma 13.21 we can view the flow fi as the sum of simple paths that 
carry flow from pi to Pi+i and simple cycles of flow 0. We computed the flow fi on a graph 
that contains G and additional edges between members of P'. If we restrict fi to the original 
graph G, then we can view fi as the sum of simple paths, each carries flow from one vertex of 
P' to another vertex of P', and of simple cycles of flow. Therefore, the restriction of fi to G is 
a flow function from a subset of P' to another subset of P' . 

Before adding fi to /, there was no residual path from P' U S to T by Invariant 15.11 From 
Lemma [32] we get that the same is true after adding fi to /. Similarly, before adding fi to / 
there was not residual path from S to P' , again from Lemma 13.21 we get that the same is true 
also after adding fi to /. We conclude that / satisfies Invariant 15.11 also after adding fi. □ 

After we added fi to /, it is possible that excess{pi) = 0, if this is the case then we are done 
with Pi. Otherwise (or if i = r and we skipped the computation of fi) we return a flow of value 
excess (pi) from pi to S and to vertices of P' with negative flow excess. We return the flow using 
a process of Johnson and Venkatesan |13j . First, we make the pseudoflow / acyclic, then we 
send back from pi a flow of value excess{pi) along a reverse topological order of /. Due to the 
antisymmetry constraint, we can view this process as adding to / a flow function f- from pi to 
a set of vertices C/ C 5* U P', each with a negative flow excess in /. This takes 0(n) time using 
the algorithm of Kaplan and Nussbaum for canceling cycles of flow [15] . 

Now the flow excess of pi is 0, this finishes the processing of pi for the case where the flow 
excess pi was initially positive. After this step P' = {pj \ j > i}. 

Lemma 5.3. After returning the excess flow from pi to S, the pseudoflow f satisfies Invari- 
ant \5.1[ 

Proof. Before returning the excess flow from pi to U Q S L) P' hy adding f- to /, there was no 
residual path from pi to any other vertex of P' , we get this from the maximality of fi, since the 
arcs with capacity which we added between vertices of P' allow to extend any residual path 
from Pi to any other vertex of P' to a residual path from pi to Pi+i. Also before we added to 
/, there was no residual path from S to P' from Invariant 15.11 Therefore, after adding /| to / 
there is still no residual path from S to P' by Lemma 13.21 (for the set A in the statement of the 
lemma we used here the set {pi} U S, note that P' does not contain pi anymore). In addition, 
before we added to /, there was no residual path from {pi} U P' U S to T from Invariant 15. H 
and this remains true after adding f- to f by Lemma 13.21 □ 

The third case of processing pi is when excess (pi) < 0. In this case we symmetrically add 
arcs {pi+i,pi+2), {Pi+2,Pi+3), ■ ■ ■ , {Pr~i,Pr), each with capacity Ai. We find a maximum flow fi, 
whose value is bounded by — excess (pi), from Pi-\.i to pi in the residual network, we remove the 
arcs we added from pi+i to pr from the graph, and let / = / + fi. Then, if excess{pi) remains 
negative, we send back flow of value —excess{pi) from T to pi. This case also keeps Invariant 
15. H the proof is symmetric to the previous case. 

When we are done, the flow excess in each vertex of P is 0, so / is a flow function. By 
Invariant 15.11 there is no residual path with respect to / from a vertex of S* to a vertex of T. 
Therefore / is a maximum flow. 

It takes us 0(n'^/^ logn) time to find fx and /y, and for each vertex of P we spend 0{n) 
time to eliminate any non-zero flow excess it has. Therefore, the total time for the procedure 
in this section is 0(n^/^ logn + rn). 



8 



We note that the bottleneck of our algorithm is the multiple-source single-sink maximum 
flow computation. For the interesting case where G is a grid, it is possible to use techniques 
similar to these presented here to improve the running time of the multiple-source single-sink 
algorithms [21 [H] to 0(n3/^). This will also improve the time bound of our multiple-source 
multiple-sink maximum flow algorithm to Oir?!"^ logn). 
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